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Abstract Weak gravitational lensing is rapidly becoming one of the principal probes 
of dark matter and dark energy in the universe. In this brief review we outline how weak 
lensing helps determine the structure of dark matter halos, measure the expansion rate 
of the universe, and distinguish between modified gravity and dark energy explanations 
for the acceleration of the universe. We also discuss requirements on the control of 
systematic errors so that the systematics do not appreciably degrade the power of 
weak lensing as a cosmological probe. 
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1 Introduction 

Within the past decade, weak gravitational lensing — slight distortions of galaxy images 
due to the bending of the light from distant galaxies by the intervening large-scale 
structure — has become one of the principal probes of cosmology. The weak lensing 
regime corresponds to the intervening surface density of matter being much smaller 
than some critical value. While weak lensing around individual massive halos was 
measured in the 1990s [BE], weak lensing by large-scale structure was eagerly expected, 
its signal predicted by theorists around the same time 3,4,5,6,7,8,9 . In this latter 
regime, the observed galaxies are slightly distorted (roughly at the 1% level) and one 
needs a large sample of foreground galaxies in order to separate the lensing effect from 
the noise represented by random orientations of galaxies. A watershed moment came 
in the year 2000 when four research groups nearly simultaneously announced the first 
detection of weak lensing by large-scale structure jlOHl 1 ifTSll 13| . Since that time, weak 
lensing has grown into an increasingly accurate and powerful probe of dark matter and 
dark energy. 
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The principal power of weak lensing comes from the fact that it responds only to 
dark matter, and not to visible (or, more generally, baryonic) matter like most other 
probes of the large-scale structure. Therefore, modeling of the visible-to-dark matter 
bias, a thorny and complicated subject, is altogether avoided when using weak lensing. 
Simulations of dark matter clustering are becoming increasingly accurate, and in prin- 
ciple there is no reason why simulation-based predictions cannot reach the accuracy 
required to model the weak lensing signal so that modeling errors do not appreciably 
contribute to the total error budget. Because of our ability to model its signal accu- 
rately, weak lensing has great intrinsic power to probe dark matter and dark energy in 
the universe. 

The other principal reason why weak lensing is powerful comes from the nature of its 
observable quantity. Galaxy shear is measured relatively straightforwardly, and comes 
from millions of galaxies typically observed in current surveys. While individual galaxy 
shear measurements do not provide cosmological information, correlation of galaxy 
shear across the sky does, especially if redshift information of galaxies is utilized to 
obtain 3-dimensional "tomography". In fact, shear-shear correlations can be used to 
obtain stringent constraints on cosmological parameters and use information from the 
survey to control the systematic errors at the same time. 

In this article, we briefly review the physics of how weak lensing probes dark matter 
and dark energy. Our emphasis is on the latter, as dark energy has emerged as the 
principal mystery of physics and cosmology (for a review of dark energy, see [14]). The 
discovery of dark energy and observations of weak lensing are both about a decade 
old, and there have been numerous investigations in recent years as to how best to use 
weak lensing to probe dark energy. This is a condensed review of the subject; for a 
more extensive treatment of weak lensing, we recommend excellent reviews by Mellier 
[15j . Bartelmann and Schneider [TB], Refregier [17| . Munshi et al. [IS], and Hoekstra 
and Jain |19j . 

2 Weak lensing basics 

2.1 Theoretical foundations 

In Newtonian Gauge the perturbed Friedmann-Robertson- Walker metric reads 

ds 2 = - (1 + 2*0 dt 2 + a 2 (t) (1 - 2$) \d X 2 + r 2 (d6 2 + sin 2 6>d<£ 2 )] (1) 

where \ i s the radial coordinate, 9 and 4> are angular coordinates, a is the scale factor, 
and we have set c = 1. Here r(x) is the comoving distance; in the special case of a fiat 
universe, r(x) = X- Finally ^ and ^ are the gravitational potentials, and $ ~ ^ in 
General Relativity and in the absence of anisotropic stresses. 

Gravitational lensing produces distortions of images of background galaxies. These 
distortions can be described as mapping between the source plane (5*) and the image 
plane (7) [20] 

Sxf = Aijdxj, (2) 
where 5x are the displacement vectors in the two planes and A is the distortion matrix 



A= ( 1 - ac - 7i -72 \ 
V ~72 1 - « + 7l / 



(3) 
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The deformation is described by the convergence k and complex shear 7 = 71 + *72 ■ 
We are interested in the weak lensing limit, where [7] <C 1. The convergence in any 
particular direction on the sky n is given by the integral along the line-of-sight 

K(fk, X )= [ X W(x')S(x')d X ', (4) 
Jo 

where 5 is the relative perturbation in matter density and 

W( X ) = ^n M H 2 g( X )(l + z) (5) 

is a weight function, iljvf is the matter energy density relative to critical, Hq is the 
Hubble constant, and z is the redshift. Furthermore 

9W =, m /;^^-^m, ( .) 

where n(x) is the radial distribution of source galaxies (normalized so that J dx n {x) = 
1) and the equality that follows the arrow holds only if all sources are at a single redshift 
z s . The distribution of source galaxies often follows a bell-shaped curve whose peak 
redshift depends on the depth of the survey. 

The quantity that is most easily determined from observations is shear, which is 
directly related to the ellipticity of the observed galaxy (in the weak lensing limit, shear 
is approximately equal to the average ellipticity). Shear is given by 

7 = 71 + ry 2 = - (ip t n - '0,22) + #,12 (7) 

where ip is the projected Newtonian potential, ip — —(1/2) J g(x) O^ij + @,ij)dx, 
and commas denote derivatives with respect to directions perpendicular to the line of 
sight. Unfortunately, this quantity is not easily related to the distribution of matter 
in the universe and the cosmological parameters. Convergence, on the other hand, is 
given by 

«= g (VUl+^,22) (8) 

which can be directly related to the distribution of matter (see Eq. (O), and is con- 
venient for comparison with theory. However, it is very difficult to measure the con- 
vergence itself, as convergence depends on the magnification of galaxies which would 
need to be measured. 

Measurements of ellipticities of distant galaxies can be used to reconstruct the "re- 
duced shear" 7/(1 — k), which is approximately equal to shear 7 in the weak lensing 
limit (for the potential deleterious effects of assuming this approximation when esti- 
mating cosmological parameters, see |21|). Measurements of shear are complicated by 
atmospheric effects such as blurring by the atmosphere, the fact that shapes of galaxies 
are irregular and not perfectly elliptical, shearing due to the instrument, etc. Correction 
of these effects and minimization of the underlying biases are very important. Recently, 
impressive attempts to find agreement between different shear estimators have been 
made in the weak lensing community with the Shear TEsting Programme (STEP 22, 
[23]). 



Fig. 1 Cosmic shear field (white ticks) superimposed on the projected mass distribution from 
a cosmological N-body simulation: overdense regions are bright, underdense regions are dark. 
Note how the shear field is correlated with the foreground mass distribution. Figure courtesy 
of T. Hamana. 

2.2 Shear/convergence power spectrum and tomography 

While individual galaxy shear (or convergence) cannot be predicted theoretically, sta- 
tistical correlations of galaxy shears can. Transforming the convergence into multipole 
space 



and assuming statistical isotropy, the power spectrum of convergence Pf is then defined 
as the harmonic transform of the two-point correlation function of the convergence 



The convergence power spectrum is identical to the shear power spectrum in the 
limit of weak distortions; P 1 {£) ~ P K (l). Either quantity is the observable quantity 
that contains a lot of information about the weak lensing field. In the limit of a Gaussian 
field, the two point function would contain all information; however, since the weak 
lensing field is nongaussian on small scales, higher-order correlations are important; 
see the next subsection. 

Weak lensing tomography — slicing of the shear signal in redshift bins - enables 
extraction of additional information from the weak lensing shear, as it makes use of the 
radial information [24] , Consider correlating shears in some redshift bin i to those in 
the redshift bin j. The tomographic cross-power spectrum for these two redshift bins, 
at a given multipole £, is defined by 




(9) 



(Ke m Kt> m <) = 5 ei e 2 8 mi m 2 P K (i). 



(10) 



k K im,i K t 'm' ,j) — &l\l-2 Sm 1 m 2 Pij(@)- 



(11) 



and can be related to theory via [4*1171124] 




(12) 



Solid: w=-1.0 
Dashed: w=-0.9 




Multipole 1 

Fig. 2 Cosmic shear angular power spectrum and statistical errors expected for LSST for 
ui = —1 and —0.9. For illustration, results are shown for source galaxies in two broad redshift 
bins, z a = 0-1 (first bin) and z s = 1-3 (second bin); the cross-power spectrum between the two 
bins (cross term) is shown without the statistical errors. We employ the redshift distribution 
of galaxies of the form n(z) oc z 2 exp(—z/zo) that peaks at 2zq = 1.0. Adopted from Ref. |14| . 



where r(z) is the comoving angular diameter distance and H{z) is the Hubble pa- 
rameter. The weights Wi are given by Wi(x) = § Om Hq Qi(x) (1 + z) where f/i(x) = 
r (x) ' dx s n l {xs)r{xs - x)/ r (Xs), and m is the comoving density of galaxies if Xs 
falls in the distance range bounded by the ith redshift bin and zero otherwise. The ob- 
served convergence auto-correlation power spectrum has additional contribution from 
the shot noise given by random galaxy shapes 

cm=pm+6 ij ^ , (13) 



where (7 int ) ' is the rms intrinsic shear in each component (which is typically around 
0.2), and fii is the average number of galaxies with well- measured shapes in the ith 
redshift bin per steradian. Clearly, the more galaxies there are and the fewer their 
intrinsic ellipticities are, the smaller the shot noise will be. 

The statistical uncertainty in measuring the shear power spectrum is [4|[7] 



AP§{t) = 



(2£+l)/ sky 



<7i 2 „t> 



(14) 



where / s k y is the fraction of sky area covered by the survey. The first term in brackets, 
which dominates on large scales, comes from cosmic variance of the mass distribution, 
and the second, shot-noise term results from both the variance in galaxy ellipticities 
( "shape noise" ) and from shape-measurement errors due to noise in the images. Fig. [2] 
shows the dependence on the dark energy of the shear power spectrum and an indication 
of the statistical errors expected for a survey such as Large Synoptic Survey Telescope 
(LSST [25]), assuming a survey area of 15,000 sq. deg. 
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2.3 Three-point correlation function of shear/convergence 

The two-point correlation function has been by far the most analyzed, and so far 
best measured, weak lensing statistic. However, higher-order correlation functions also 
contain significant cosmological information; chief among these is the three-point cor- 
relation function. With purely Gaussian shear or convergence field, the three-point 
function is of course zero (within measurement and cosmic- variance errors). While 
primordial nongaussianity can produce a nonzero three-point function, nonlinear clus- 
tering of structure presents the most widely anticipated — and in fact guaranteed — 
signal in the three-point function 6 26 ,27 ,28 ,29 30 . 

The three-point correlation function of the convergence in multipole space can be 
related to the bispectrum of the convergence via 

(^rrnNfemj^ms) = ^ ^ j B W 3 ( 15 ) 

where the term in parentheses is the Wigner 3j symbol. The convergence bispectrum 
is essentially a projection along the line of sight of the matter overdensity bispectrum 

RK / (2li + l)(2l 2 + l) (2*3 + 1) / h l 2 h \ 

B ^ - V 4^ ) 

where W(z) is the same weight function as in Eq. (|12p . The (convergence) bispectrum is 
defined only if the following relations are satisfied: \£j — £k\ < £% < \£j+£k\ for{i, J, fc} £ 
{1,2,3} and £\ + £2 + £3 is even. To compute the bispectrum of the convergence, 
therefore, we need to supply the matter bispectrum B(k\, k2, k$, z) which, on mildly 
nonlinear scales of interest, can be evaluated using the halo model or perturbation 
theory techniques, or else calibrated directly from cosmological N-body simulations. 

Measurements (and detections) of the three-point function have been first made in 
the VIRMOS-DESCART survey [311132], and, more recently, in the COSMOS survey 
[33] . Much better constraints are anticipated from future wide-field surveys |34] , How- 
ever, including the covariance between the power spectrum and the bispectrum — a 
difficult quantity to estimate in its own right — may degrade the combined constraints 
relative to a naive (i.e. uncorrelated) combination of their constraints |35j . 



3 Galaxy-galaxy lensing and galaxy cluster counts 

One effective application of weak lensing is to measure the correlation of shear of the 
background galaxies with mass of the foreground galaxies. This method goes under the 
name of "galaxy-galaxy lensing" 2 , 36 , 37 , 38 , 39 , and essentially measures the galaxy- 
shear correlation function across the sky. Galaxy-galaxy lensing measures the surface 
mass density contrast A£(R) 

AE(R) = E{< R) - S(R) = r crit x lt (R). (17) 

S(< R) is the mean surface density within proper radius R, £(R) is the azimuthally 
averaged surface density at radius R (e.g. [3,40 ), jt is the tangentially projected shear, 
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Fig. 3 Surface density profile measurements obtained from galaxy groups in the COSMOS 
survey |44II45| . The legend shows various components of the halo model fit to the profile. Figure 
kindly provided by Alexie Leauthaud. 



and the critical surface density S CT n is a known function of the distances to the source 
and the lens. The surface density can be further inverted to obtain radial density 
profiles of dark matter halos [41,42 . More accurate modeling of the galaxy-galaxy 
lensing signal, that takes into account clustering with other halos and the baryonic 
content of each halo, can be done using the halo model of large scale structure [43| . 
Figure [3] shows the surface density profile of a sample of ~600 galaxies in the COSMOS 
survey with stellar masses of order ~ 1O 11 M0 and at z ~ 0.5, together with halo model 
fits [44U45| . Current measurements constrain the density profiles and bias of dark matter 
halos 41,46,47,48 , and the relation between their masses and luminosities [491150] . In 
the future, galaxy-shear correlations have potential to constrain dark energy models 
[51j and modified gravity models for the accelerating universe [52] . 

Weak lensing signal can also be used to detect and count massive halos — partic- 
ularly galaxy clusters. This method, pioneered recently [531154] , can be used to obtain 
cluster samples whose masses are reliably determined, avoiding the arguably more dif- 
ficult signal-to- mass conversions required with the X-ray or optical observations [55 , 56 
57 58]. Much important information about the dark matter and gas content of galaxy 
clusters can be inferred with the combined lensing, X-ray and optical observations. 
This has recently been demonstrating with observations of the "Bullet" cluster [59] , 
where the dark matter distribution inferred from weak lensing is clearly offset from the 
hot gas inferred from the X-ray observations (see Fig. [4]), indicating the presence and 
distinctive fingerprints of dark matter. 

Finally, one can use the abundance clusters detected through weak lensing to con- 
strain cosmological parameters (e.g. [601161] ). The number of halos of in the mass range 
[M, M + dM] in a patch of the sky with solid angle dfl and in a redshift interval 
[z, z + dz] is given by 

d 2 N r\z)dn{M,z) 

dQdz-W) dM dM (18) 

where r(z) and H(z) are the comoving distance and the expansion rate respectively, 
and n(M,z) is the number density (the "mass function") that can be calibrated from 
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Fig. 4 X-ray emission from the "Bullet" cluster of galaxies observed by ground and space 
telescopes. Colored features correspond to X-ray emission observed by Chandra space telescope, 
while the green contours correspond to mass reconstruction from weak lensing observations. 
The "bullet" is a smaller subcluster which has passed through the larger cluster and whose 
hot gas is seen in X-rays, being clearly distinct from the total mass of the system. Adopted 
from Ref. [59]. 

numerical simulations. Because weak lensing signal is directly sensitive to the cluster 
mass, the usual uncertainty in mapping from the observable signal (X-ray, or optical 
light, etc) and mass is largely avoided. However, contamination of the cluster shear by 
the projected mass (i.e. large-scale structures between us and the cluster) can seriously 
degrade the effectiveness of lensing-detected cluster counts (e.g. 62 ). It is therefore nec- 
essary to use N-body simulations to calibrate purity (contribution of false detections) 
and completeness (fraction of detections relative to the truth) of the shear peaks 63,64, 
65,66,67,68,69 , and the jury is still out as to whether the projection systematics can 
be controlled to sufficient accuracy. An alternative to using the weak-lensing-detected 
clusters is to apply weak lensing mass measurements to clusters detected via X-ray [70] 
or optical [7T] observations. 



4 Cosmological constraints 

First detections of weak lensing by large-scale structure have been announced nearly a 
decade ago [1011111 112 13 . Their results were in mutual agreement and consistent with 
theoretical expectations, which is remarkable given that they were obtained indepen- 
dently. Meanwhile, additional measurements have been collected using ground-based 
72 , 73 , 74 , 75 , 76 , 77 78 and space-based 79 , 80 , 44 , 81 observations, leading to increas- 
ingly better understanding of systematic errors. 



4.1 Current constraints 

Current ground-based weak lensing surveys cover of order ~ 100 square degrees. Space- 
based surveys currently cover only a few square degrees, but utilize significantly more 
distant galaxies whose lensing signal is more pronounced. 
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Fig. 5 Example of cosmological constraints from weak lensing data, adopted from Canada- 
France-Hawaii Telescope Legacy Survey (CFHTLS; 84 ). Left panel: E and B mode mea- 
surements; note that B modes are consistent with zero, and their scatter is consistent with 
expectations from the alignment between intrinsic shapes of galaxies. Right panel: Constraints 
(68% and 95% confidence) in the Qm-cs plane. Blue contours are from weak lensing, green con- 
tours are from the CMB (WMAP 3-year; 85 ), and the orange contours are the two combined. 
Adopted from Ref. [HI]. 



Modest sky coverage of current surveys, combined with evolving understanding of 
the systematic errors, means that only a few parameter combinations can be measured 
to an interesting accuracy as of this writing. In particular, the combination agil^f, 
where erg is the amplitude of mass fluctuations, is measured to about 5-10% from 
the ground [821178] and ~15% from space |83| . and is in concordance with values of 
Qm ~ 0.25 and erg ~ 0.8. An example adopted from Ref. [84] is shown in Fig. [5] 
with measurements of E and B mode correlation functions of shear in real space (left 
panel) , and the resulting constraints on ag and Q^i (right panel) . Such constraints are 
impressive given that first detections of weak lensing by large-scale structure have been 
made less than a decade ago. 

Weak lensing also probes dark energy, though at the present time it needs to be 
combined with other cosmological probes in order to produce interesting constraints 
on the equation of state w. Recent work illustrates how weak lensing helps tighten 
the constraints from SNe la and CMB [86] mostly due to its sensitivity to growth of 
density perturbations. It is also possible to place constraints on specific dark energy 
model from current weak lensing data [87]. Upcoming weak lensing surveys from the 
ground and space hold promise to measure dark energy significantly better, as we now 
discuss. 



4.2 Future constraints 

The errors in the cosmological constraints can then be easily and straightforwardly 
forecasted using the Fisher matrix formalism. Consider a set of cosmological parameters 
Pi (i = 1, . . . , iVpar) and weak lensing observations that result in measurements of the 
shear power spectrum Cj k (£), defined in Eq. (|13[) and where j and k run over the 
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Fig. 6 Illustration of forecast constraints on dark energy parameters, adopted from Ref. I14| . 
Shown are 68% C.L. uncertainties for one version of the proposed SNAP experiment, which 
combines a narrow-area survey of 2000 SNe to z = 1.7 and a weak lcnsing survey of 1000 sq. 
deg. Left panel: Constraints in the Om- w plane, assuming constant w; the vertical axis can 
also be interpreted as the pivot value w p for a time-varying equation of state. Right panel: 
Constraints in the wo-w a plane for time-varying dark energy equation of state, marginalized 
over I2]vi for a flat Universe. 



tomographic redshift bins. Then the Fisher matrix is defined as 




where C is the column matrix of the observed power spectra and Cov 1 is the inverse 
of the covariance matrix between the observed power spectra whose elements are given 
by 

Cov [C%(£), Cfi(0] = (2 f+T)/ sky [C&WM + ■ (20) 

assuming that the convergence field is Gaussian, which is a good approximation at 
I < 3000. By the Kramer-Rao inequality, cosmological parameters can typically be mea- 
sured to no better than cr(pj) = 1/y/Fa (unmarginalized error), or er(pj) = \J {F~ 1 )a 
(marginalized over all other parameters). In practice, these idealistic error bars are usu- 
ally a good approximation to the truth if the parameters are sufficiently well measured 
so that their covariances are ellipsoid-shaped and not "banana-shaped" regions. 

Weak lensing measurements can determine cosmological parameters to an excellent 
accuracy due to the fact that modeling the weak lensing signal only requires modeling 
of dark matter, though baryonic collapse at centers of dark matter halos and galaxy 
alignments complicate the picture somewhat and contribute to systematic errors dis- 
cussed below. 

First cosmological parameter accuracy projections from weak lensing have been 
made in [88] while increasingly more sophisticated estimates 89,90,34,91,92,93,94, 
95 96 indicate prospects for determining the parameter erg to much better than 1% 
and the equation of state of dark energy to a few percent accuracy in the future. An 
example of such projections is shown in Fig. [SJ where we show constraints on Qm 
and constant w (left panel) and wq and w a (right panel) from SuperNova/Acceleration 
probe (SNAP 97 ), a proposed space-based experiment to probe dark energy. These 
constraints marginalize over ~ 5 other cosmological parameters, and allow for the 
presence of systematic as well as statistical errors. 
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Fig. 7 Example of how weak lensing can constrain modifications of gravity parametrized by 
the growth index 7, adopted from Ref. 11031 . Left panel: Auto-correlation power spectra of 4-bin 
weak lensing tomography for the two values of the growth index, with statistical errors shown 
for the 7 = 0.55 model. This shows that a relatively small signal due to modified gravity 
can be detected using future weak lensing observations. Right panel: Bias in constraints in 
cosmological parameters wo and w a if changes in modified gravity were ignored, for the = 
0.1 modification relative to the fiducial model. Recall, the growth index in DGP braneworld 
gravity differs from the General Relativity value by Aj ss 0.67 — 0.55 = 0.12. 

Finally, we mention the proposed technique to correlate the mass of foreground 
galaxies with the shear of the background galaxy population, and thus isolate geometric 
distance factors in photometric redshift bins. This approach sometimes goes under the 
name of "cross-correlation cosmography" [98], and is very closely related to galaxy- 
galaxy lensing. Cross-correlation cosmography is potentially less sensitive to systematic 
errors as well as assumptions about the nonlinear clustering of dark matter than the 
shear-shear correlation measurements, and can significantly complement cosmological 
constraints obtained from the latter. [98U99II100U 101U 102U51] 



5 Weak Lensing and Modified Gravity 

Weak lensing is a particularly effective discriminator of modified gravity explanations 
for the accelerating universe. Modified gravity stipulates that the reason for the ap- 
parent acceleration are modifications to General Relativity on large scales (e.g. in 
the infrared). General Relativity has to be preserved on the solar-system and galactic 
scales in order to satisfy stringent observational constraints, but modifications on large 
scales could masquerade as an apparent acceleration when the equations of standard 
Einstein's gravity are used to interpret the observations. 

Gravitational lensing is particularly useful probe of modified gravity because it 
probes the sum of the two gravitational potentials <f> + & (see Eq. |T])), while particle 
dynamics probes & alone (e.g. [10411105] !. Since modifications to gravity typically affect 
the potentials differently, combination of weak lensing with other cosmological probes 
can in principle be used to differentiate modified gravity from dark energy. For a 
detailed review of how weak lensing can be used to test gravity, see the review |106] in 
this Special Issue. 

The simplest approach to distinguish between modified gravity and honest-to-god 
dark energy is to parametrize anomalous growth of density perturbations; the most 
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popular approach is via the "growth index" parameter |107| . The linear growth factor 
g(a) = 5 (a) /a) (scaling out the matter dominated universe behavior 5 oc a), where a 
is the scale factor, is approximated by 

g (a) = ef° dhla[aM(a) "'- :L] , (21) 

where 7 is a new parameter called the growth index (not to be confused with shear!). 
This fitting function is accurate to 0.3% compared to the exact solution within general 
relativity for a wide variety of physical dark energy equation of state ratios including 
equations of state w(a) = wq + w a (l — a), provided that the function takes the value 
7 = 0.55 + 0.05(1 + w(z = 1)) [TOT]. The fitting formula also fits [108] the linear 
growth in the popular "braneworld gravity" model of Dvali, Gabadadze and Porrati 
(DGP 109 ), but with a different value of the growth index, 7 ~ 0.67. More generally, 
measuring 7 can distinguish between modified gravity and dark energy even if the 
two predict identical expansion history (e.g. distances in the universe; for a similar 
approach not using the growth index, see [110II111] !. Future expectations are for 7 to 
be measured to accuracy of about 0.02-0.03 from weak lensing data combined with 
other cosmological probes that measure the distance scale (SNe la, CMB) or growth 
(cluster abundance) 103, 112 ; sec Fig. [7] for an example. 

Weak lensing can also probe specific models where dark matter or dark energy 
interactions are modified. For example, models of massive gravity, where the Poisson 
equation is phenomenologically modified with the Yukawa term, have been constrained 
with weak lensing observations [11311114] ; preliminary constraints have been imposed 
on some other models as well [115] . In the future, weak lensing data will enable a much 
more ambitious battery of tests of modified gravity, essentially enabling comparison 
of distance and growth functions at multiple redshifts and on multiple angular scales 



6 Second-order effects as a signal for cosmology 

So far we have discussed weak lensing of galaxies by the large-scale structure, a principal 
signature of weak lensing used in cosmology. However, objects other than galaxies are 
also weakly lensed, and their properties can sometimes be used to our advantage. 

In particular, high-z type la supernovae (SN la) will be gravitationally lensed by 
the foreground large-scale structure. The useful property of SN la is that they have 
nearly uniform intrinsic luminosities. Weak lensing will magnify SN la and perturb 
these luminosities, and the resulting magnifications will be correlated and may be 
measurable with future surveys. The resulting magnification of supernova luminosity 
distance is 

H=[(l- - H 2 ]" 1 ~ 1 + 2k + 3k 2 + | 7 | 2 + ... , (22) 

where k is the lensing convergence and I7I = \/ "ff + To is the total shear at the 
angular location of the SN. Each supernova's luminosity L and distance will be 
perturbed as SL/L — 2\Sd^/d^\ = fi— 1, where jj, is magnification at that location and 
out to redshift of the SN. The magnification of SN — really, the correlated pattern 
of the departures of their observed distances from the true values — comes for free 
with standard measurements of SN distances to measure the expansion history of the 
universe. 
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Measurements of cosmic magnification of supernovae therefore complement galaxy 
shear measurements in providing a direct measure of clustering of the dark matter [120| . 
As the number of supernovae is typically much smaller than the number of galaxies, the 
two-point correlation function of lensed SN la is more noisy than, and not as effective 
as, the galaxy shear. However, weak lensing of SN la probes magnification directly, and 
can be used to measure the second-order corrections in Eq. (|22|l between magnification 
and convergence (or equivalently, shear and convergence) 120 . Additionally, the power 
spectrum of magnification can be used to measure the cosmological parameters such 
as as |121] , While this application of weak lensing is rather challenging and requires a 
large number of SN la to be effective, we emphasize that it comes for free in wide-field 
surveys that measure galaxy shear as well as SN la, and combines weak lensing shear 
and magnification information in the same field, providing a number of cross checks on 
the systematics. 

7 Systematic errors and their calibration requirements 

Because weak lensing has so much statistical power, it is imperative to control the 
systematic errors so that they do not appreciably degrade the statistical errors on the 
weak lensing power spectrum (see e.g. the small error bars in Fig. [5]), and the similarly 
small errors on the cosmological parameters (see Fig. [6]). 

7.1 Observational errors 

One of the most important systematic effects on weak lensing are the photometric 
redshift errors — uncertainties in measuring redshifts of sheared galaxies. For weak 
lensing measurements that simply integrate signal along the line of sight, redshift errors 
are clearly not critically important, as only the overall radial distribution of galaxies, 
n(z), needs to be known (see Eq. [6]). However, for tomographic measurements, the 
photometric redshift errors lead to leakage of galaxies from their original redshift bin 
and contamination of signal in the bin where they are misplaced. 

The required accuracy of photometric redshift bins has been studied in detail |122| ; 
and it was found that the mean and spread of the relation between photometric and 
(true) spectroscopic redshift, averaged over bins of Az = 0.1, both need to be known to 
better than about 0.003 |122| : see the right panel of Fig. [8] While errors in individual 
galaxy redshift are greater than this by about an order of magnitude, it is reasonable 
to expect that modern photometric redshift techniques (e.g. |123j l can achieve this 
required accuracy in each redshift bin. Moreover, a possibility of catastrophic redshift 
errors — cases where the true redshift is misestimated by a significant amount — is 
equally troubling, and requires knowledge of outlier "islands" in the 2spec--z p hot plane 
to a similar sub-percent accuracy |1241I125] . 

Large size of current and future surveys, with millions and up to a billion of ob- 
served galaxies, requires the use of photometric redshifts. These photometric redshifts 
will be calibrated by obtaining spectra of a subset of galaxies, and "training" the 
photometric spectra using the spectroscopic subset. Of order 10 5 spectra seem to be 
required for future surveys [122,93,126,124 , although this number strongly depends 
on the expected scatter in the photometric redshifts around the true value and other 
details. Moreover, the spectroscopic galaxies need to be a representative sample of the 
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Fig. 8 Requirements on the accuracy of photometric redshifts to perform weak lensing tomog- 
raphy, adopted from Rcf. 122 . Contours show error degradations in constant w (left panel) 
and those in w a (when where both wq and w a are varied; right panel), as a function of prior 
knowledge of the mean (x-axis) and scatter (y-axis) in the relation between the true (spectro- 
scopic) and photometric redshift. The degradations are defined as actual errors in cosmological 
parameters relative to the errors with perfectly known photometric redshifts. For example, the 
right panel shows that no more than 50% degradation (factor 1.5 on the contours) in the error 
in w a requires knowledge of the photometric redshift mean and scatter to better than ~ 0.003 
in each redshift bin. 



whole population — for example, very high redshift galaxies in the survey need to be 
represented in the spectroscopic subsample — and this adds another layer of challenge 
to the program of photometric redshift calibration. 

Errors on the observed shapes of galaxies are manifold, and include imperfect knowl- 
edge of the point-spread function, atmospheric blurring, thermal and mechanical effects 
of the telescope, etc. More generally, errors in measurements of galaxy shear can be 
described as either additive or multiplicative errors |127j . The measured shear of source 
galaxies at redshift z s and in direction n, 7(z s ,n), can be written as 

7(z s , n) = 7 (z s , n) x [1 + fi(z s , n)] + 7 add (z s , n) (23) 

where 7(z s ,n) is the true shear, fi(z s ,n) is the multiplicative and 7 a dd(- z si n ) the 
additive error ( 127 ; see also 128 ). Multiplicative errors can be caused by, for example, 
finite size of the point-spread function with which the galaxy image is convolved while 
additive errors can be created, for example, by atmospheric blurring. The requirements 
are stringent; for example, future surveys need to calibrate the mean multiplicative 
errors to about 0.1% (i.e. Sfi w 0.001) accuracy in each of ~ 10 redshift bins in 
order that dark energy constraints not be significantly degraded |127] . However, what 
really helps weak lensing is the possibility of "self-calibrating" the systematic errors 
— determining a reasonable set of the systematic error nuisance parameters from the 
survey concurrently with the cosmological parameters without appreciable degradation 
in accuracy on the latter [127 129, 130 . With self-calibration, the survey itself is used 
to partially calibrate the systematic effects. 
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7.2 Theoretical modeling systematics 

Predicting the weak lensing shear power spectrum requires knowledge of the dark mat- 
ter power spectrum (see Eq. ll2[) ; intrinsic uncertainty in predicting the power spectrum 
in the nonlinear regime (on scales I > 200, or below about 1 degree) constitutes the 
principal theoretical weak lensing systematic. Ref. [131) found that power spectrum 
on scales of ~ 1 Mpc (corresponding roughly to a few arcminutes on the sky), where 
weak lensing measurements are most accurate, needs to be calibrated to a few percent. 
While not trivial, this accuracy seems within reach due to recent impressive progress 
in obtaining precision results from large numerical simulations (e.g. |132] V 

One additional difficulty, currently only partially explored, is how the addition of 
gas to pure dark matter simulations affects the results; it appears that baryonic physics 
can be phenomenologically modeled with sufficient accuracy 133,134,135 . Another 
challenge to theoretical predictions is that a large suite of numerical simulations appears 
necessary to cover the multi-dimensional cosmological parameter space; fortunately, 
clever statistical techniques can be employed to cover the parameter space so that a 
very reasonable total number of simulations (~ 100) is sufficient |136II137] . 

Weak lensing provides a large amount of information, given by shapes and redshifts 
of millions (or, in the future, billions) of galaxies. Therefore, there is redundancy in 
weak lensing measurements that can be used to minimize the effect of systematic 
errors with negligible degradation in the statistical errors on cosmological parameters. 
For example, by taking only the cross-power spectra for cosmology (PQ (I) with i ^ j 
in Eq. JT2}), one can minimize biases due to correlations between intrinsic galaxy 
ellipticities [138] (see also [139] ). Similarly, one can perform "nulling tomography" of 
weak lensing to selectively throw out the small-scale information in the convergence 
power spectrum that is most sensitive to the unwanted biases [140] , 

8 Conclusions 

Weak lensing probes dark matter and energy in the universe in multiple ways. With 
galaxy-galaxy lensing, it probes the structure of individual dark matter halos. With 
detecting galaxy clusters from the shear of the foreground galaxies, it enables the 
cosmological number-count test that is sensitive to distances and growth of density 
perturbations, and thus to temporal evolution of dark energy. And with measuring 
the correlation of shear signal across the sky, weak lensing again probes distances and 
growth, but now with a completely different set of statistical and systematic errors. 

In this review article, we briefly covered the aforementioned cosmological uses of 
weak lensing; we did not mention some related topics that are covered elsewhere in this 
Special Issue (for example, lensing of the CMB |141| ). We also emphasized the intrinsic 
power of weak lensing, and the potential of weak lensing to self-calibrate an impressive 
range of systematic errors and still produce excellent cosmological constraints. This 
situation can perhaps be contrasted with strong gravitational lensing, where only a 
modest number of multiply imaged systems (~ 100) is currently available, and where 
it is important to have very precise a-priori constraints on the systematic effects in 
order to extract the cosmologically interesting information. 

Current efforts in weak lensing research emphasize performing observations with 
ever wider and deeper surveys, improving the accuracy and reliability of weak lensing 
measurements, improving the understanding of nonlinear clustering of dark matter 
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and baryons with N-body simulations, and devising clever methods to self-calibrate 
the systematic errors and measure new properties of dark energy and modified gravity. 
If these efforts continue at the current rate of progress, weak lensing is likely to firmly 
establish itself as one of the most reliable and powerful probes of dark energy and 
matter in the universe. 
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